Controllable Raman-like nonlinearities from non-stationary cascaded quadratic 

processes 
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We show that useful non-instantaneous nonlinear phase shifts can be obtained 
from cascaded quadratic processes in the presence of group velocity mismatch. 
The two-field nature of the process permits responses that can be effectively 
advanced or retarded in time with respect to one of the fields. There is an 
analogy to a generalized Raman-scattering effect, permitting both red and blue 
shifts of short pulses. We expect this capability to have many applications in 
short-pulse generation and propagation, such as the compensation of Raman- 
induced effects and high-quality pulse compression, which we discuss. 
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1. Introduction 



In the last decade, there has been much interest in the 
nonlinear phase shifts produced by the cascaded inter- 
actions of two or three waves in quadratic (x' 2 '') nonlin- 
ear media. Large nonlinear phase shifts of controllable 
sign can be generated, and numerous applications can 
be envisioned on the basis of such a capability 1 . The 
prototypical quadratic process is second-harmonic gener- 
ation (SHG). During the propagation of a fundamental- 
frequency (FF) field along with its second-harmonic 
(SH), the FF accumulates a nonlinear phase shift 
if the process is not phase-matched. With long pulses 
(nanosecond-duration in practice) the FF and SH fields 
overlap temporally despite their different group veloc- 
ities. In this so-called stationary limit, an effective 
Kerr nonlinearity is obtained (except at high inten- 
sity, when the fundamental field is depleted), and this 
can be a surrogate for the bound-electronic cubic (x^) 
nonlinearity^. The cascade nonlinear phase shift can be 
thought of as arising from an effective nonlinear refrac- 
tive index, i. e. the real part of an effective susceptibility. 
The residual SHG that occurs in the phase-mismatched 
process can similarly be considered the analog of two- 
photon absorption (the corresponding imaginary part of 
the effective susceptibility). 

The use of cascaded quadratic processes with ultra- 
short pulses is complicated substantially by the group- 
velocity mismatch (GVM) between the FF and SH 
fields^. GVM causes the fields to move apart in time, 
which reduces their coupling, and thus the magnitude 
of the cascade effects. In addition, the temporal profile 
of the nonlinear phase shift becomes distorted. Devia- 
tions of <fr NL (t) from the pulse intensity profile hamper 
or preclude applications that involve soliton-like pulse 
shaping. The solution to this problem amounts to re- 
covery of the stationary regime: For given value of the 
GVM, the phase mismatch is increased so that the cycles 
of conversion and back-conversion that generate the non- 



linear phase shift occur before the pulses move apart from 
each other in time. Liu et al. showed that acceptable 
phase-shift quality can be obtained if at least 2 conver- 
sion cycles occur per characteristic GVM length Lqvm = 
0.6cTo/(n gi i — n g _2), which implies Afc > 47t/Lgvm^- 
Here, c is the speed of light in vacuum, tq is the full- 
width at half-maximum (FWHM) of the pulse, and n Sj i 
and n g> 2 are the group refractive indices for the FF and 
SH, respectively. In the limit of large phase mismatch an 
exact replica of a cubic nonlinearity is asymptotically ob- 
tained^. The disadvantage of working with large phase 
mismatch is reduced magnitude of the nonlinear phase 
shift. GVM thus places a strong constraint on applica- 
tions of cascade phase shifts. As one example, increas- 
ing GVM reduces the fraction of launched pulse energy 
that evolves into a soliton, eventually to zerot To date, 
applications of cascade phase shifts with femtosecond 
pyjgggEjjSjgjio ]j ave a n been demonstrated under station- 
ary conditions. Approaching the stationary boundary, 
GVM coupled with self-phase modulation has been ob- 
served to asymmetrically broaden the pulse spectrur 
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The nonlinear refraction experienced by an ultrashort 
pulse in a cubic nonlinear medium arises predominantly 
from bound-electronic and nuclear (i.e. Raman) contri- 
butions to the nonlinear response. Here we show that 
cascade phase shifts produced under non-stationary con- 
ditions mimic the Raman response, with some remark- 
able properties. Frequency shifts of controllable sign and 
magnitude can be impressed on short pulses. These ef- 
fective Stokes and anti-Stokes Raman processes complete 
the analogy between cascade nonlinearities and true cu- 
bic nonlinearities, while maintaining the new degree of 
freedom provided by the quadratic interaction - control 
of the process through the phase mismatch. An interest- 
ing feature of the non-stationary cascade process is that 
it provides a controllable non-instantaneous (and there- 
fore nonlocal) nonlinearity. That the GVM alters the 
quadratic processes and produces deviations from a Kerr 
nonlinearity is well-known^. However, to date these ef- 
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fects have been perceived as distortions to be avoided. 
Just as the ability to control the sign and magnitude 
of an effective nonlinear index has enabled a new class 
of applications^, controllable Raman-like processes can 
be expected to create substantial new opportunities. We 
can think of cascade nonlinear processes with short pulses 
as dividing naturally into two classes, separated by the 
criterion for obtaining nonlinear phase shifts that mimic 
those of purely electronic origin. With this view, half of 
the possibilities have yet to be explored. Some examples 
will be discussed. 

2. Analytical Approach 

The Kerr-like nonlinearity that arises from the x : x' 2 ' 
process in the stationary limit can be understood qualita- 
tively as follows. A small fraction of the FF is converted 
to the SH, which accumulates a phase difference before it 
is back-converted. The nonlinear phase shift impressed 
on the FF is delayed by one full cycle of conversion and 
back-conversion. However, as long as the SH is not dis- 
placed temporally from the FF, the phase shift on the 
FF will be proportional to its intensity profile (Fig. QJ. 
With short pulses, the GVM becomes important if the 
fields separate by approximately the pulse duration be- 
fore a cycle of conversion and back-conversion is com- 
plete. (We assume that the effect of the difference be- 
tween the group-velocity dispersions (GVD) of the FF 
and the SH is negligible; in practice, it typically is much 
weaker than the inter-pulse GVM. This will be discussed 
quantitatively below.) Thus, after one cycle of conversion 
and back-conversion, the intensity profile of the SH field 
is retarded or advanced with respect to that of the FF, 
depending on the sign of the GVM. The corresponding 
delay of the nonlinear effect is slightly smaller or larger 
than one full conversion cycle. As a result, an effectively 
advanced or retarded phase shift is accumulated by the 
FF. The corresponding effect in the spectral domain is 
a frequency shift toward the blue or the red. Such fre- 
quency shifts have been predicted through numerical cal- 
culations^. 





Fig. 1. An illustration of the cascaded quadratic pro- 
cesses under phase-mismatched conditions. The FF is 
partially converted to the SH and then back-converted. 
Dashed (solid) lines are for the case of zero (nonzero) 
GVM. 

The propagation of the FF and the SH are governed 



by coupled equations within the slowly varying envelope 
approximation (SVEA)— . We neglect self- and cross- 
phase modulation due to x'' 3 '' , consider only the tempo- 
ral dimension, and assume conditions for type I second 
harmonic generation, but the results can easily be gener- 
alized. 
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Here ai and a 2 are the normalized FF and SH field 
amplitudes. Time is normalized to the initial pulse du- 
ration, t = t/ro, and the scaled propagation coordi- 
nate is £ = z/Lqvm- Here, 5j = Lqvm/Lds.j, where 
Lds.j — 0.322to 2 /GVD(cJj) are the dispersion lengths 
with j = 1,2. The parameter (3 — AfcLcvM where 
Afc = k 2uJ — 2k w is the normalized FF-SH wave-vector 
mismatch. 

Consider the simple but common case when only FF 
light is incident on the quadratic medium. In the limit 
of large phase-mismatch, conversion and back-conversion 
occur rapidly and most of the energy resides in the FF 
at all times. A relation between the FF and the SH 
amplitudes can be derived as an expansion in powers of 
j3~. By eliminating a 2 in Eq. [21 and keeping terms up to 
order (1//3 3 ), an equation of motion for the FF field can 
be derived: 
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The first three terms constitute a nonlinear 
Schrodinger equation (NLSE), which is generalized 
by the fourth and the fifth terms that describe the 
effects of the GVM, and the mismatch of the GVD of 
the FF and SH, respectively. The final term is negligible 
(at order 1//3 2 ) since it is proportional to the square of 
the first derivative of the field envelope (a small quantity 
in the SVEA). In addition it is even in time for well 
behaved fields (ai(t)), and hence cannot contribute to 
the frequency shifting process. We neglect the GVD 
mismatch since its effect is much smaller than inter-pulse 
GVM for typical nonlinear media under the assumed 
conditions of small conversion to the SH. In that case, 
Eq. reduces to: 
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We notice that Eq. ^resembles the Chen-Liu-Lee equa- 
tion (CLLE), which is integrablsi^. The difference is the 
presence of a cubic nonlinear term. It can be shown that 
Eq. El reduces to the CLLE: 



.dq{Z,T) d 2 q{Z,T) 



dZ 



dT 2 



= 0, (5) 



with the substitution, 



T ) = co exp[i(ci£ + c 2 t)]q(£, r), 



(6) 



where |c | 2 = B 2 ^j%J%, Ci = (Si/3 2 /8, and c 2 = (3/2, 
followed by the coordinate transformation 



z = t 



(7) 



Hence, Eq. 0] is integrable as well. We note that the 
CLLE is further related to the well-known, integrable 
derivative nonlinear Schrodinger equation, via a Gauge- 
invariant transformation^. 

To further understand the effects of the lowest order 
correction from GVM in Eq. we can compare it to 
the equation governing the propagation of a single field 
envelope a\ under similar approximations to Eas. 11121 but 
for a Kerr-nonlinear material with finite Raman-response 
time T K *2£ 



0. (8) 




NLSE 



Here 7 = U20Jq / (ckw 2 ) for a Gaussian beam with of 
frequency ujq and waist w. Comparing Eq. [3] to Eq. 0] 
we see that they are similar but with the Raman term of 
Eq. El replaced by Tf |a x \ 2 da x /dT where Tf = -2i//3. 
While the correspondence is not exact since the func- 
tional dependence is different (~ a\d\a\\ 2 / dr for Raman- 
scattering vs. ~ \a\\ 2 da\/dT for the cascaded process), 
some qualitative understanding can be gained from con- 
sidering the effective cascaded response with T£. ff ~ i/fi: 
first, the cascaded correction is imaginary and hence does 
not contribute directly to the phase (unlike the Raman 
response). Rather it alters the field envelope. The en- 
velope change subsequently couples to the phase profile 
through the remaining terms of Eq. ^ so the frequency 
shift occurs through a higher order process. Second, 
saturates with 1/(3 ~ l/Afc(A), unlike the Raman re- 
sponse which does not depend strongly on wavelength. 



This saturation will be explored in greater detail in Sec- 
tion [3J Note that the effective response for the cascaded 
process can be > two orders of magnitude greater than 
that of Raman-scattering, so that significant frequency 
shifting is possible in cm of quadratic material (vs. me- 
ters of fiber with Raman). 

For a qualitative understanding of the effect of the 
GVM term, we decompose the field a%(£, r) in Eq.QJinto 
its amplitude and phase with the substitution di(£, r) = 
t) exp(i0(£, t)), where m(£,t) and 4>(£,t) are real 
functions. The evolution of the amplitude and the phase 
is then given by 



du du d<f> Si d 2 <p 

di = 1 d^dr + Y u d^ 



2 2 du 



(9) 
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For the sake of simplicity, we concentrate on the non- 
linear terms and ignore the dispersion, which amounts to 
neglecting the terms with higher-order time derivatives. 
With this simplification, we obtain 



du 



,du 



3 



2 U Q T . 
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(10) 

(11) 



If we assume (3 is large, t) and </>(£, r) can be cal- 
culated by expanding in powers of 4, similar to the pro- 
cedure used to obtain Eq. Keeping terms up to order 
jp, Eg. ITUI shows that r) = u(t) + O(-jp), that the 
field amplitude is approximately unchanged. Thus, inte- 
gration of Ea.lTTl yields 



1 ~, 1 



«^, T ) = -~^( r )£ + 0(_). 



(12) 



Substitution of this relation back into Eq. ^] gives 



9 4 o du 

1 -(P U d-T- 



(13) 



The first term on the right represents the Kerr-like 
nonlinear phase shift and the second term corresponds 
to the non-instantaneous nonlinear response due to large 
GVM. 

The relation in Eq. 1131 provides a valid description of 
the phase evolution only in its early stages, before the 
field amplitude is modified significantly, and in the ab- 
sence of dispersion. Within these approximations, the 
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effect of the GVM on the nonlinear phase shift can be 
illustrated for a given pulse shape: ai(0,r) = sech(r). 
Integration of Ea. 1 131 yields 



1 , 2 2 



sech^(r)(l + i— sech(r) tanh(r))£ + 0o,(14) 
P 



where </>o is an integration constant. The temporal 
asymmetry of the GVM contribution shifts the peak of 
the nonlinear phase shift. 

Likewise, Fourier transforming to the frequency do- 
main, the contributions of the Kerr-like and the GVM 
terms in Eq. 0] can be calculated for the pulse shape to 
be 



\a 1 \ 2 a 1 + 2i^-\a 1 \ 2(> "' 



1 



Ot 







-(1 + oj 2 )(3 + 2cj//3)sech(7rw/2), 



(15) 



which has a bipolar shape. Positive frequency compo- 
nents are attenuated and negative frequencies are ampli- 
fied, or vice versa, depending on the sign of the phase- 
mismatch to GVM ratio. This result is expected to hold 
in general for any smooth, single-peaked pulse shape for 
which ai(0,r) — > for |t| — > oo. Such a frequency shift 
is expected from the GVM term, which has an odd-order 
time derivative. 

3. Numerical Analysis 

Although the approximate one-field equation (Eq. [3} is 
useful for a qualitative understanding, it is necessary to 
consider the coupled equations (Eqs. I1I2|) for a quanti- 
tative description. To this end, we numerically solve a 
version of Eas. 11121 that has been generalized to include 
self- and cross-phase modulation terms due to the cubic 
nonlinearity. We use a different field normalization here 
to facilitate comparison to experimental parameters. 

.dAi Zi d 2 A 1 

2 ■ — 

dz 2Lds,i 9t 2 
+ A* 1 A 2 exp(iAk(Z I z)) 



'NLA 



|A 1 | 2 + 2|^1 2 | 2 )A 1 =0, 



(16) 
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dz 2Lds,2 dr 2 Lgvm Ot 
] -A\ exp(-iAfc(Z/z)) 



n(cj%) , 



n{u> 2 ) 
n{uji) Z 



n(u 2 ) L NL 



I —(2\A 1 \ 2 + \A 2 \ 2 )A 2 =Q. (17) 



Here the FF and SH envelopes (Ai and A 2 , respec- 
tively) are in units of the initial peak FF field Aq (related 



to the initial peak FF intensity by Iq = 



72) 



and n 2 is the Kerr nonlinear index. The cubic nonlinear 
length characterizing the pulse propagation is -kjvij = 
c/u>jn 2 I (the length over which the accumulated nonlin- 
ear phase shift is 1) for frequency loj with j = 1,2. The 
propagation coordinate z is normalized to the quadratic 
nonlinear length Z\ — nAi/27rd gAo which characterizes 
the strength of the nonlinear coupling. d c ff is the effec- 
tive quadratic nonlinear coefficient. Time r, L£>s,j 1 and 
£gvm are defined as for Eas. 11121 



Eqs. 11611 71 are solved using a symmetric split-step 
beam propagation method^. The simulations assume 
typical conditions for femtosecond pulses in quadratic 
nonlinear crystals. As an example, we calculate the 
propagation of 120-fs pulses with a peak intensity of 
50 GW/cm 2 in a 10 cm-long quadratic medium. The 
launched pulse shape is chosen to be Gaussian and the 
center wavelength is 790 nm. The quadratic medium used 
is barium metaborate (Ba 2 B04 or BBO), for which the 
material parameters are: n = 1.63, d e g — 1.82 pm/V, 
n 2 = 3.2 x 10~ 16 cm 2 /W, GVM = -186 fs/mm, the 
FF (SH) GVD = 70 (190) fs 2 /mm, and the FF (SH) 
third-order dispersion is 50 (81) fs 3 /mm. Note that the 
true cubic nonlinearity (n 2 ) is included in the calcula- 
tions. Two-photon absorption is neglected since it is 
small for BBO at 790 nm. The phase-mismatch is set 
to be Afc = 57r/mm, corresponding to a self-defocusing 
nonlinearity and a magnitude that is about half the min- 
imum value to obtain a Kerr-like phase shift. This par- 
ticular set of conditions, except for the crystal length, is 
chosen to correspond to experiments that are described 
later in the text. As expected, the spectrum of the pulse 
shifts to higher frequencies as it propagates through the 
quadratic medium. The evolution of the spectrum is 
shown in Fig.^a). Initially, the frequency-shift increases 
linearly with propagation distance, but eventually the 
process saturates (Fig. Etb) ) . The spectrum of the SH 
field shifts opposite from the FF (i.e. to lower frequen- 
cies), prior to saturation. This saturation is expected, 
since the effective response of the cascaded process de- 
creases (or is distorted) with increases frequency shift. 
For this choice of pulse parameters and phase-mismatch, 
the saturation begins beyond 3 cm, which is close to the 
maximum length of available BBO crystals. The cubic 
electronic nonlinearity of the quadratic material is in- 
cluded here for complete correspondence with experimen- 
tal parameters, however the close agreement in Fig. |2{b) 
between the saturation trend with and without n 2 in- 
dicates the dominance of the quadratic process in the 
frequency shifting dynamics. In Fig[2^b), the presence 
of cubic nonlinearity slightly reduces the resulting fre- 
quency shift, as expected due to its self-focusing phase. 
One might expect the material's Raman response to be 
relevant to the frequency shifting dynamics studied here, 
however the Raman response of BBO with ^100- fs pulses 
is small compared to the cascaded response. In compari- 
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son to Raman, the cubic electronic nonlinearity included 
above in Fig. |^b) is a larger effect, even though it al- 
ters the frequency shifting process indirectly through the 
nonlinear phase. 

Frequency shift (FWHM) 



-2 2 4 6 




z (mm) 

Fig. 2. (a) Evolution of the spectrum along the propa- 
gation direction. Shift is in units of the initial spectral 
FWHM, ~ 3.7 THz. The scale bar shows spectral inten- 
sity in arb. units, (b) Weighted average frequency shift 
as a function of propagation distance. Dashed line indi- 
cates fit to region of linear shift. Dash-dots show similar 
results in the absence of \^ ( n 2 = 0). 

The pulse propagation is dominated by an interplay 
of GVD and the effective nonlinearity from the cascaded 
process in the form of soliton-like dynamics. In the time 
domain, the pulse undergoes compression since the en- 
ergy is more than the amount needed to balance the dis- 
persive effects. The intensity profiles before the onset of 
saturation and at the exit face of the crystal are plotted in 
Fig.[3]along with those of the launched pulse. The FF un- 
dergoes a steady compression accompanied by energy loss 
to the SH: At z = 12 mm, its FWHM is 110 fs with 64% 
of the pulse energy remaining in the FF. With propaga- 
tion, the pulse shape becomes slightly asymmetrical. The 
asymmetrical structure develops into a secondary pulse 
in the final stages of the propagation, which corresponds 
to the secondary structure of the spectrum (Fig. |2Ja)). 
At z — 100 mm, the FWHM of the main peak is reduced 
to 40 fs while ~ 36% of the launched energy is retained 
in the FF. The temporal profiles are displaced since the 
pulse experiences different group velocities as its central 
frequency changes. 

A similar picture emerges for the total frequency 
shift for fixed propagation distance and varying phase- 
mismatch. The magnitude of the nonlinear phase, and 
hence the frequency shift, is proportional to 1/A.k be- 
fore it saturates. Loss to SH conversion increases with 
decreasing |Afc|, so there exists a trade-off between the 
magnitude of the frequency shift and loss. We define a 
figure-of-merit (FOM) for the shifting process as the ra- 
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Fig. 3. Intensity profiles of the FF at z = mm (dash- 
dotted line), z = 12 mm (solid line), z — 100 mm 
(dashed-line). For the launched pulse Lds.i = 74 mm. 

tio of frequency-shift to energy content in the SH field, 
which attains a maximum for phase-mismatch values 
slightly below that of the minimum for obtaining Kerr- 
like nonlinear phase shifts (Fig. ■ This is demonstrated 
in Fig. which shows simulations of 100-fs, 200-pJ 
pulses with center wavelength 1550 nm (and peak in- 
tensity ~ 5 GW/cm 2 ). The material parameters used 
correspond to the quadratic material periodically-poled 
lithium niobate (PPLN): n = 2.14, d cS = 16.5 pm/V, 
ri 2 = 3.2 x 10" 15 cm 2 /W, GVM = -370 fs/mm, and 
the FF (SH) GVD = 100 (400) fs 2 /mm. Under these 
conditions, the stationary boundary for Kerr- like phase 
shifts corresponds to |Afc| > 257r/mm. Notice that much 
larger frequency shifts can be generated closer to phase- 
matching, but with larger SH-conversion. 

The non-instantaneous nature of the cascaded 
quadratic process with significant walk-off between the 
FF and the SH is demonstrated by the nonlinear phase 
shift imposed on the FF. Simulations confirm the afore- 
mentioned expectations: Effectively retarded or ad- 
vanced phase shifts are imposed on the FF depending 
on the sign of the GVM (Fig.|5j). 

4. Experimental Observation of the Frequency- 
shift 

Experiments were performed with 120-fs, 0.6-mJ pulses, 
centered at 790 nm, generated by a Ti:sapphire regener- 
ative amplifier. The launched pulse shape was approxi- 
mately Gaussian with a clean spatial profile and the peak 
intensity is estimated to be 50 GW/cm 2 . A 17 mm- 
long piece of BBO served as the quadratic medium. The 
GVM length is Lqvm = 0.38 mm, for which the crite- 
rion for a Kerr-like phase shift implies Afc > 10.47r/mm. 
Both blue and red shifts are experimentally available 
via positive and negative phase-mismatch, respectively. 
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Fig. 4. Frequency shift (crosses) and FOM (circles) as 
a function of phase-mismatch. Similarly to Fig. [21 fre- 
quency shift is measured in units of the initial FWHM 
(here ~ 4.4 THz). Note that GVD is chosen to be normal 
(anomalous) for Afc > (Afc < 0) to support soliton-likc 
pulses. 




Time (FWHM) 

Fig. 5. The phase impressed on the FF for zero (solid 
line), positive (dashed line), and negative (dash-dotted 
line) GVM. 



temporal profile of the pulse did not change significantly 
in these experiments. The experimental results are com- 
pared to the results of numerical simulations which con- 
tained no free parameters and were based on experimen- 
tal conditions. In Fig. the calculated spectrum of the 
unshifted pulse (367r/mm) is normalized so that it con- 
tains the same power as the measured unshifted spectrum 
and all other traces have the same relative scaling, so that 
the units of all the given spectra are the same. With this 
in mind, there is excellent agreement between the mea- 
sured spectra and the simulations. In particular, the shift 
increases greatly from Afc = 197r/mm to Ak = 57r/mm. 
The latter is the phase-mismatch for which the ratio of 
the spectral-shift to SH conversion should peak. The de- 
pendence of the frequency-shift on the phase-mismatch, 
as summarized by the inset of Fig. El agrees qualitatively 
with the results of Fig. however, the phase-mismatch 
corresponding to maximum spectral-shift to SH conver- 
sion is different than in Fig. 0] as a consequence of differ- 
ent physical parameters in the experiment. 





10 20 30 40 
Ak(:t mm" 1 ) 



760 770 780 790 800 810 820 830 
Wavelength (nm) 

Fig. 6. Experimental (solid lines) and simulated 
(dashed lines) spectra for phase-mismatches of 57r/mm 
and 367r/mm. The latter serves as control. Inset: Exper- 
imental (symbols) and calculated (solid line) frequency 
shift for different values of phase-mismatch. As in Fig. El 
frequency shift is measured in units of the initial FWHM, 
- 3.7 THz. 



However, red-shifts occur with self-focusing nonlinearity, 
which limits the peak intensity available without contin- 
uum generation and crystal damage. Consequently we 
focus here on blue-shifts. 

Increasing frequency-shift of the FF was observed with 
decreasing phase- mismatch (Fig. [HJ . The inset of Fig. [Bl 
shows the spectral shift versus phase-mismatch for self- 
defocussing phase shifts, and the main figure shows exam- 
ple spectra. The data presented for Afc = 367r/mm serve 
as a control experiment: At such a large phase mismatch 
the cascade nonlinear phase and the Kerr nonlinearity 
are negligible and the spectrum is indistinguishable from 
the spectrum of the launched pulses (not shown). The 



5. Applications 

We have shown that cascaded quadratic processes 
under phase-mismatched conditions and in the pres- 
ence of significant GVM (typical conditions with 
femtosecond-duration pulses) result in an effectively non- 
instantaneous cubic nonlinearity. The response time 
is controllable by appropriate choice of the phase- 
mismatch. In addition to providing an intuitive picture 
for the effect of phase-mismatch on the propagation of 
femtosecond pulses, this nonlinear process offers some 
unique features. The nonlocal nature of the Raman-like 
cascade nonlinearity is interesting in its own right. Non- 
locality of the cubic nonlinearity has been shown to ar- 
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rest self-focusing collapse and to stabilize solitons, for 
exampl e - 17 ! 18 . The nonlocal nature of the cascade pro- 
cess under non-stationary conditions can be controlled 
or tailored to specific situations via the phase-mismatch. 

Many applications of a controllable effective Raman 
process can be envisioned. Perhaps the most obvious one 
is the cancellation of the Raman shift that a short pulse 
accumulates as it propagates in optical fiber. For exam- 
ple, in telecommunication systems with bit rates above 
~ 20 Gbit/s, the pulse duration is short enough that 
timing jitter is dominated by jitter arising from Raman- 
induced frequency shifts^. 

In high-energy short-pulse fiber amplifiers, the nonlin- 
ear phase shift can be controlled reasonably well by the 
technique of chirped-pulse amplification, and as a result 
an equally important limitation to pulse energy is stim- 
ulated Raman scattering— The red-shifts produced by 
Raman scattering can be compensated by blue-shifting 
the pulses prior to, or following, propagation in fiber. 
As an example, we calculate the pre-compensation of the 
Raman-induced red-shift of a 100-fs, transform-limited 
pulse centered at 1550 nm in standard single-mode fiber 
(modal area of 80 fim 2 and GVD of -23 ps 2 /km). The 
pulse energy is 1 nJ. The quadratic medium is a 4 cm- 
long wave-guide written in PPLN. The modal area of 
the waveguide is 40 /im 2 and the GVD for the FF is 
100 ps 2 /km2i. The phase-mismatch is set to Afc = 
207r/mm. The pulse is first blue-shifted in the PPLN 
waveguide and then propagates in the fiber. These cal- 
culations indicate that the central wavelength can be kept 
at 1550 nm, following propagation in up to 50 cm of fiber. 
If no pre-compensation is utilized, the pulse is red-shifted 
to 1800 nm (Fig. 0. This result nicely complements the 
previously-established conclusion that the cascade non- 
linearity can be used to compensate the nonlinear phase 
shift produced by the electronic Kerr nonlinearity under 
similar conditions— 

Other potential applications include devices that con- 
vert peak power to frequency-shift, which can be used to 
switch wavelength channels or intensity discrimination 
with the addition of a frequency filter—. 

We consider the application to pulse compression in 
some detail. For pulse energies in excess of 1 mJ, meth- 
ods based on cubic nonlinearity for the generation of 
extra bandwidth cannot be used, owing to the limita- 
tions of excessive nonlinearity in single-mode waveguides 
and material damage through self-focusing for unguided 
geometries. Self-defocussing nonlinearities in quadratic 
media address these difficultie a 5 i 24 . The generalization 
of this approach to include frequency shifts as described 
here enables us to implement an analog of Raman-soliton 
compression 2 ^: high-order solitons are formed, producing 
a compressed primary pulse that undergoes a continuous 
self- frequency shift. An advantage of this approach is 
that the pedestal commonly produced by Raman-soliton 
compression consists mainly of unshifted frequency com- 
ponents. These components can be eliminated with a 
frequency filter to yield a pedestal-free pulse. 




Wavelength (\im) 

Fig. 7. Pulse spectrum after propagation in fiber with- 
out pre-compensation (dash-dotted line) and after cas- 
cade pre-compensation stage (dashed line) and subse- 
quent propagation through fiber (solid line). Dots in- 
dicate the launched pulse spectrum. 




200 

Time Delay (fs) 

Fig. 8. Temporal profile of compressed pulses before 
(dashed line) and after spectral filtering (solid line) of 
the unshifted frequencies. Inset displays the compressed 
pulse spectrum before (dashed line) and after filtering 
(solid line) . Dash-dotted lines indicate the launched tem- 
poral profile/spectrum. 



Numerical simulations for realistic parameters demon- 
strate the utility of this approach: 100-fs, 0.6-nJ pulses 
with Afc = 507r/mm compress to 20-fs upon propaga- 
tion through a 6 cm- long waveguide in PPLN (Fig. [HJ. 
The pulse quality (Q c ), defined as the ratio of energy 
contained within the FWHM of the pulse to that of the 
initial pulse, is calculated to be 0.65. The unshifted com- 
ponents can be filtered out to produce a longer (38 fs) 
but much cleaner (Q c = 0.91) pulse. Compression in a 
second 2.5-cm-long PPLN crystal generates a 15 fs pulse 
with virtually no additional degradation in pulse qual- 
ity (Qc — 0.90). The resulting pulse after two stages of 
compression contains ^50% of the launched pulse energy. 
Similarly, calculations indicate that compression factors 
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of up to 3 should be attainable with 1 mJ pulses and us- 
ing a bulk BBO crystal at 800 nm, and experiments are 
underway to verify this compression. In addition to the 
high pulse quality, a practical advantage of this approach 
is that larger nonlinear phase shifts can be produced at 
the smaller phase mismatches needed in comparison to 
compression in quadratic media under nearly-stationary 
conditions. 

6. Conclusion 

In summary we have demonstrated a new capability of 
cascaded quadratic processes under phase-mismatched 
conditions: Effectively retarded or advanced nonlinear 
phase shifts can be impressed on a pulse in the presence 
of significant GVM between the FF and SH frequen- 
cies. The frequency-domain manifestation of this non- 
instantaneous nonlinear response is red- or blue-shifts of 
the pulse spectrum. The direction and the magnitude 
of the frequency-shift is controllable by the choice of the 
phase-mismatch. Just as effectively instantaneous phase 
shifts from cascaded processes are analogous to bounded 
electronic (x*' 3 '') nonlinearities for negligible GVM, these 
non-instantaneous phase shifts in the presence of strong 
GVM are analogous to nuclear (Raman-induced) nonlin- 
earities. 

We expect the unique features of these processes to 
find many applications. Here, we numerically demon- 
strated compensation of Raman-induced frequency-shifts 
and high-quality pulse compression assuming typical con- 
ditions for femtosecond pulses in common quadratic non- 
linear media. More generally, however, spectral shifts 
from cascaded quadratic processes should be applica- 
ble to all processes involving Raman-induced frequency- 
shifts, but with the added freedom of sign and magnitude 
control. 
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